Strong-coupling surface polaron in a magnetic field 



Hua Xiukun 1 



Department of Physics 
Suzhou University, Suzhou, 215006, China 



Abstract 

By applying variational method, the strong electron-surface optical(SO) phonon interaction 
and the weak electron-longitudinal optical (LO) phonon interaction are studied systematically 
in present paper. The formula of induced energy V e - so , V e -i Q in both the ground and excited 
states are given. The numerical results show that as the distance between the electron and 
the polar crystal surface increases, the electron-LO phonon interaction energy increases yet the 
electron-SO phonon interaction energy decreases. The numerical results also show that both 
electron-LO and electron-SO phonon interaction are enhanced as the magnetic field strength 
increases. 
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1. Introduction 

The behavior of polaron has received much attention in the past decades [1,2]. Huybrechts 
[3] proposed a linear combination operator method and investigated the property of polaron in 
strong electron-phonon coupling case. Tokuda [4] added another variational parameter to the 
momentum operator and studied the ground-state energy and effective mass of the bulk polaron. 
Gu [5] discussed the the ground-state energy, the first-excited state energy and the effective 
mass of intermediate-coupling polaron in a polaron slab. When external magnetic field is applied 
to the electron-phonon interaction system, the property of magnetopolaron attracted many 
investigators' interests. Whitfield et al [6] studied the effectiveness of adiabatic approximation 
theory for a polaron in a magnetic field and they found that the usual adiabatic theory applied 
to a polaron in a strong magnetic field does not give the right weak coupling limit. Ercelebi 
[7] studied the two-dimensional magnetopolaron in the strong-coupling regime and proposed a 
modified coherent phonon state. Hai, Peeters et al [8] did research into the polaron-cyclotron- 
resonance spectrum resulting from interface- and slab-phonon modes in a GaAs/AlAs quantum 
well and an experiment was done to investigate the resonant magnetopolaron effects due to 
interface phonons in GaAs/AlGaAs multiple quantum well structures by Wang, Nickle et al 
[9]. Recently, the author etc have studied the electron self-energy and effective mass by the 
powerful Green's function method [10]. 

For the bulk polaron, the perturbation theory is applicable when on < 1 and the LLP (Lee- 
Low-Pines) variational method is effective when an < 6. For the surface polaron, Pan [11] 
showed that when the electron-surface optical(SO) phonon coupling constant a s > 2.5, the 
strong coupling theory must be applied. Therefor ,in some polar crystals, the weak electron-LO 
phonon coupling and strong electron-SO phonon coupling should be considered together. So 
far, the research work in such a field has been scarce 

In the present paper, using variational method and taking weak electron-LO phonon and 
strong electron-SO phonon interaction into account, we give the ground and excited states 
energies of the system, the electron-LO phonon interaction induced energy V e -ho an d the 
electron-SO phonon interaction induced energy V e so- Using polar crystal AgBr as an example, 
numerical evaluation are also presented. 

2. Theory 

Considering the following system: the semi-infinite space z > is occupied by the AgBr 
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crystal, whereas the space z < is a vacuum. The external magnetic field is applied along the 
z direction and an electron is moving in the polar crystal. B=(0,0,B), choosing the symmetric 
gauge for the vector potential, i.e. A=B/2(-y,x,0), the Hamiltonian describing an electron 
coupled to both LO phonons and SO phonons can be written &s(h = m b =l, m b is the band 
mass of the electron) 

H = H + H± + H e - ph = H± + Hj i , (1) 
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In the above equations, a\{ak) is the creation(annihilation) operator of bulk LO phonons 
with three-dimensional wave vector k, b\{bt) is the corresponding operator for the SO phonons 
with two-dimensional wave vector q. P/j — (p x ,p y ) and p = (x,y) are the electron momentum 
and position vector in xy plane, respectively, and p z is the electron momentum in z direction. 
k// = (k x , k y ) is phonon wave vector in the xy plane. lo\ and cu so are the frequencies of the bulk 
LO and SO phonons. eo and are the static and high-frequencies dielectric constant of the 
crystal respectively. A and V are the surface area and the volume of the crystal. 

The polaron system wave function can be separated into the electron and the phonon parts. 

|V) = l0e}|<M (7) 

With \4> P h) = W1W2IO), |0) is the phonon ground state, and the explicit form of U\,U2 will be 
introduced in the following. 
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Let us start from the unitary transformation: 

ui = exp(-i(J2 a^lakk// ■ P + J2 a 2&jM • p)) (8) 
k q 

With ai = 1 and 02 = corresponding to the weak-coupling limit case and the strong-coupling 
limit case respectively. The parts of polaron Hamiltonian can be transformed into: 

H> = u^H oUl = 1 -{P U - £ a\a k k n f + ^(x 2 + y 2 ) + ^L z 

k 

+-7rJ2 a l a k(ykx ~ xk y ) + ^ui a\a k + ^u so b\b q , (9) 

1 k k q 

H' e _ lo = u^ l H e _ loUl = J2[V k *sin(zk z )al + h.c.} , (10) 

k 

K-so = u^He-soii! = J2[V* exp(-qz)exp(-iq ■ p)b\ + h.c] , (11) 

q 

H' ± = u^H ±Ul = H ± . (12) 

Following the scheme of Huybrechts, we introduce the creation and annihilation operators 
cj and Cj by 

2 

i 

Where the subscript j refers to the x and y directions, A is the variational parameter, and c](c 3 -) 
is Boson operator. 

Rewriting (9), (11) using (13), (14), one gets: 

H o = ljE( c H + c i c i) + + c i c i)] ~^r^ a l a kk//j(c] + C j) 

3 3 k,j 

1 1 UJ 2 

+ 7; S a l a kk 2 / + - alalaka^kf/k'// + -| ^(c,ct + c ) Cj ) 

j k q 

+ ^^J2 a k a k[K(c y -cl) -k y (c x -cl)} , (15) 

#e- so = Y:[V q *exp(-qz)exp(Y: ^(c,- - cj))&t + , (16) 



Pi = — ( c i + c i) ' ( 13 ) 
Pi = 4f(9 " cj) . (14) 



Yet, H' e _ lo and H' ± are invariable in forms. Let continue to do the u 2 transformation: 

u 2 = exp[J2(a{f k e-^rPo _ a k f* k e ik //^) + "£(b\g q e^ - b q g* q e*<*)] (17) 
k q 

Where fk(fk) an d 9q(g q ) are variational parameters, and po = (xo, Uo) is the electron orbit center 

where x = x/2 — p y /cu c , y Q = y/2 + p x /u c .The necessity of making the phonon deformation 

centered at p was emphasized in an elaborate discussion by Whitfield et al [6]. 

u^a k u 2 = a k + f k e^ k /rP\ u 2 l a\u 2 = a\ + f k e ik 'r(», (18) 
u 2 \u 2 = b q + g q e-^ ) u 2 x b\u 2 = b\ + g q e^ . (19) 

Performing u 2 transformation and using (18) and (19), one can get Hq, H"_ lo and H"_ so 
after lengthy calculation. 
(l)The ground state: 

Determining the ground state |0 > in the new representation by 

Cj \0 >= a k \0 >= b q \0 >= (20) 
one gets the energy E'j, (0) of the ground state 

e'^o) = <hj,>=<h» + hi Io + h:-„> 

= \ + f x + E^IM 2 + Y,"so\g q \ 2 + \ E \h\ 2 k 2 „ + 1 E \h\ 2 \h>\ 2 k ■ k> 

k Q k k,k' 



k 1 A 

+ J2[V k *sin(zk z )f* k + V k sin{zk z )f k ]exp{-^-{- + — )) 
k 8 X u 2 c 

-1 



a 1 A 
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Minimizing E"/(0) with respect to f k (f k ) and g q (g* q ) and neglecting higher order of f k , one 
gets 

Ej / (0) = ± + £ + V e -U0) + V e -„(0) (22) 
Where V e _z o (0) and V e - so (0) are the correction energy induced by electron-LO phonon and 
electron-SO phonon interaction respectivelywith 

\V k \ 2 sin\zk z )exp{-^{\ + ±)) 

k 2 k //+ WO 
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\V q fexp(-2 q z)ex P (-i({ + ±)) 
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Considering the explicit forms of 14 and V q , (23) and (24) can be reduced to 
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In eq.(22) ,if neglecting V r e _; o (0),V r e _ so (0) and minimizing -Ey/(0) with respect to A, one can 
obtain A = cj c and £//(0)=|u; c . Obviously, this is the electron ground state energy of Landau 
Level. 

The effective Hamiltonian of the polaron system in the ground state H e ff is 

H eff = E'; / (0) + H ± = ^- + ^ + ^- + V e . lo + V e -so + V img , (27) 



with 
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(2) TTie excited state: 

The excited state 1 1 > in the new representation is 

|1 >= c ]\0 >= 4|0 >= b\\0 > (28) 

considering 

< lle^-i, . >= - - + g)] (29) 



similarly, one gets 



£7/(1) = ^ + M + Wjo + Wso + X4_ /o (l) + V e - so (l) (30) 



Where V e _; (l) and V e - so {\) are the correction energy induced by electron-LO phonon and 
electron-SO phonon interaction of the excited state respectivelywith 
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3. Discussion and Conclusion 

Taking the polaron in polar crystal AgBr as an example, we perform a numerical eval- 
uation. The parameters concerned are listed as follows: £ =10.6, £oo=4.68, hwi Q — 17 .lmeV ', 
hu> so — 1 A. 5meV , ai— 1.56, a s =2.56, m fe =0.22m, m is the free electron mass. All the above pa- 
rameters are taken from [12]. 

Fig.l Fig.2 

Fig.l shows the relationships between the induced potential |V e _; o (0)| resulting from the 
electron-LO phonon interaction, the induced potential |14_ so (0)| resulting from electron-SO 
phonon interaction, and their sum |T4_; o (0)| + |V e _ so (0)| and the coordinate z. From fig.l we can 
see that the induced potential V e -i o (0) increase with increasing coordinate z, whereas V e _ so (0) 
decrease with increasing coordinate z. Compared with the weak electron-phonon interaction 
case [13], | V^»_ so (0) | and |V e _i o (0)| are larger than the corresponding weak coupling ones,yet 
their lines tendency is similar to the weak coupling ones. When z < 10A, |V e _ so (0)| is larger 
than |V r e _; o (0)|, whereas for z > 35 A, |V e _ so (0)| becomes smaller compared with |V r e _i o (0)| and 
so |T4_j o (0)| is the dominant term. |V e _; o (0)| + |V e _ so (0)| decreases with increasing z. When z > 
35 A, \V e -i o (0)\ + |T4_ SO (0)| approaches |V e _; o (0)| because |T4_ SO (0)| is much smaller compared 
with \V e _ lo (0)\. 

Fig.2 show the induced potentials as a function of the electron coordinate z for different 
magnetic fields. In fig.2, the dot-line stands for the case when B=4T and the solid-line stands 
for the B=8T case. From fig.2, we can see that |l / e _i o (0)| and |V e _ so (0)| for B=8T are larger than 
the corresponding ones for B=4T. This result show that external magnetic field can enhance the 
interaction between an electron and phonons. Let us give a possible physical explanation of the 
result. When magnetic field is applied, the polarizability of the crystal increases and the density 
of phonons increases. Because the strength of the electron-phonon interaction is proportional 
to the number of phonons, the interaction energies, for both |V e _; o (0)| and |V e _ so (0)|, increase 
with increasing magnetic field strength. 

Fig.3 
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Fig. 3 illustrates the energy difference induced by the electron-LO (SO) phonon interaction 
between the excited and ground states . From fig. 3 one can see that the difference between 
V e - so (l) and V e - so (0) decreases , yet the difference between V e -i (l) and V e -i o (0) increases as 
the coordinate z increases. As z is small (z < 10A) the magnitude of V e - so {\) —V e - 8O (0) is much 
large than that of V e -i (l) — V e -i o (0), and when z > 33A \V e - so (l) — V e _ so (0)| becomes less 
than | V e _ ?0 (l) — K-2o(0)|. These results give further implication that the electron-SO phonon 
interaction play a dominant role as electron are near the polar crystal surface. 

In summary by applying variational method, the strong electron-surface optical(SO) phonon 
interaction and the weak electron-longitudinal optical (LO) phonon interaction are studied sys- 
tematically in present paper. The formula of induced energies V e _ so (0), V e _; o (0),V r e _ so (l),V e _/ o (l) 
are given. The numerical results show that as the distance between the electron and the polar 
crystal surface increases, the electron-LO phonon interaction energy increases yet the electron- 
SO phonon interaction energy decreases. The numerical results also show that both electron-LO 
and electron-SO phonon interaction are enhanced as the magnetic field strength increases. 
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Figure Captions 



Fig.l. |V e _ so (0)|, | V^_/ o (0) | , and |T4_ SO (0)| + |V e _; o (0)| as a function of coordinate z. 

Fig. 2. |V e _ so (0)| and |V e _; o (0)| as a function of coordinate z for different magnetic fields. The 
solid line and the dot line represent B=8T and B=4T, respectively. 

Fig. 3 \V e - ph (l) — V e - ph (0)\ as a function of coordinate z at B=10T. 
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